Centrality dependence of the thermal excitation-energy deposition in 8—15 GeV/c 

hadron-Au reactions 
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The excitation energy per residue nucleon (E*/A) and fast and thermal light particle multiplicities 
are studied as a function of centrality defined as the number of grey tracks emitted N SIcy and by 
the mean number of primary hadron- nucleon scatterings ({u)) and mean impact parameter ({6}) 
extracted from it. The value of E*/A and the multiplicities show an increase with centrality for 
all systems, 14.6 GeV p-Au and 8.0 GeV 7r-Au and p-Au collisions, and the excitation energy per 
residue nucleon exhibits a uniform dependence on N glay . 

PACS numbers: 25.40.a, 25.40.+t, 25.70.Pq 



I. INTRODUCTION 



Multifragmentation is the process by which a heavy nucleus excited through nuclear collision decays into nuclear 
fragments with a range of atomic masses and energies. First observed in cosmic ray emulsion experiments as a 
starburst of grey and black emulsion tracks emanating from a single incident track, it has since been studied using 
beams of projectiles including single hadrons (pions, protons, and antiprotons), light ions, and heavy ions accelerated 
■"sj" ■ to energies ranging from a few hundred MeV up to Tevatron energies of 350 GeV Q, I i. 

Analysis of the hadron-induced reactions suggests a fast evolutionary mechanism that is dominated by two stages [1, 
5]. Initially the system is heated by a prompt nuclear cascade that produces high-energy, forward peaked particles, 
commonly identified as grey particles, named according to their appearance in emulsion. The system then cools 
via preequilibrium particle emission, leading to the second major stage, an excited, thermalized nuclear residue that 
• ■ subsequently undergoes statistical decay. Analyses of the second stage emission products indicate that for events 
that deposit the highest excitation energies, the residue undergoes multifragment decay on a near-instantaneous time 
scale 0, consistent with a nuclear liquid-gas phase transition and perhaps critical behavior 0, 

The process of multifragmentation has been thoroughly investigated as a function of incident beam energy [3[ and 
projectile species but a systematic study of its centrality dependence has not yet been performed. In heavy ion col- 
lisions, the experimental definition of centrality usually refers to the mean number of primary collisions or the impact 
, parameter distribution extracted from forward particle multiplicities using a standardized set of assumptions, and it 
is often used to characterize nuclear collisions at GeV energies and above [l(| • Centrality analyses often provide 
useful insights into particle production mechanisms and enable systematic comparisons with other experiments and 
theoretical models, where comparable measures of centrality can be precisely determined. The ISiS multifragmenta- 
tion program at the Brookhaven National Laboratory Alternating Gradient Synchrotron (BNL-AGS) with its large 
acceptance and event-by-event charged particle identification is ideally suited to study the centrality dependence of 
multifragmentation. In this paper we present results from the first multifragmentation centrality analysis using the 
data from the ISiS experiments E900 and E900a. 
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II. EXPERIMENT 



The experimental grey particle and excitation energy distributions were obtained at the AGS using secondary 
beams of untagged 14.6 GeV/c protons and tagged beams of 8.0 GeV/c ir~ and antiprotons incident on a 197 Au 
target. The highest energy data sets for each beam species were selected to optimize the Glauber model assumptions 
required for the centrality analysis. Reaction products were measured with the ISiS detector array, which consists of 



162 triple detector telescopes arranged in a spherical geometry and covering 74% of 47r [ll|. Detector telescopes were 
composed of a low-pressure gas ionization chamber, followed by a 500 /im silicon detector for measuring low-energy 
fragments and a 28 mm Csl crystal for detecting more energetic particles. Telescopes provided charge identification for 
1.0<E/A<90 MeV reaction products and particle identification in the range 8.0<E/A<92 MeV, primarily hydrogen 
and helium isotopes. 

Of particular relevance to the present study, the kinetic energy of particles that punched through the Csl crystal was 
derived from the Csl energy loss in the Si-Csl particle-identification spectrum and energy-loss tables. This procedure 
permitted identification of particles up to 350 MeV in energy, which were assumed to be primarily hydrogen nuclei and 
were classified as grey particles. Full description of the experimental details, data analysis, calorimetry procedures 
and error estimates are given in 0, [l2j . 



III. ANALYSIS 

We follow the procedure developed by BNL-E910 [HI, 0, EH, m which the conditional probability for the grey 
track distribution as function of the number of hadron-nucleon scatterings (v), P(N grcy \v), was determined by a fit 
to the measured grey particle distributions. Here we use a modified exponential form for P(N gTey \i>), 



N gtcy {v) = Cl [1 - (exp(- i ,/c 2 )) C0 ] . (1) 

and assume a binomial distribution for the number of grey tracks given Z total protons. Eq.[T] preserves the dominant 
linear term of the second order polynomial used in (l3j , but accommodates the N gley saturation in a more natural way. 
7Vg re y was defined as the number of identified protons with 30<E<350 MeV. This energy range is meant to separate 
fast and thermal protons, and is consistent with previous definitions of N gTcy . The lower energy cut is identical to the 
one used by BNL-E910, although their high energy cut of 585 MeV is higher than that used by most experiments. 
The dependence on the low energy cut is included in the study of systematic errors. The full functional form is given 
byEq.[H 

P(N grcy ) = ]TP( V grc »7r(z,), (2) 

V 

\ jv grcy J 
= VgrcyQ) 

z ' 

The 7r(i/) Glauber distribution for the number of primary hadron-nucleon interactions is calculated with a Monte Carlo 
optical model [l6j using a Wood-Saxon nuclear density profile assuming free inelastic cross-sections for protons at 14.6 
GeV/c and pions and antiprotons at 8.0 GeV/c, given in the particle data book as 30, 20, and 47 mb, respectively [l7l ]. 
The Glauber distribution for the E900 data set includes a 10% admixture of 14.6 GeV/c pions [T^, also using a 20 mb 
cross-section. Fits of Eq. [2] to the fast proton N gTcy distributions are shown as solid lines in Fig. Q] panels (a), (b), 
and (c). The x 2 / n df obtained for the proton, pion, and antiproton systems were 30.2, 11.9, and 1.0 respectively. 
To estimate systematic errors, fits were also performed with a version of Eq. [2] with Co set to unity (dashed) and to 
the second order polynomial (dotted) used in [13| . The ^ jndj values for these fits were orders of magnitude larger 
for some systems. In a study of systematic errors reported in [l3[ in which the N grcy analysis was compared to a 
microscopic cascade RQMD model [l8[ the overall systematic error was estimated to be 10-20%, with the dominant 
contribution coming from the lower momentum bound in the N gvey definition. For this reason we also performed fits 
to the combined fast and thermal proton distributions 8<E<350 MeV, shown in the Fig. [T] panels (g),(h), and (i), 
and to the sum of the fast proton, deuteron, and triton distributions, shown in the Fig. [1] panels (d), (e), and (f). 
Note for the fast deuteron and triton definitions, the low energy cut is set to 49 MeV rather than 30 MeV used for 
protons. The systematic error estimate also includes fits with cross-sections varied by ±2 mb. The values of (v) and 
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FIG. 1: (Color online) Fits to grey track distributions for incident protons are shown at the top, negative pions below, and 
antiprotons at the bottom. Fits are shown for iVgrcy distributions with three definitions: fast protons (a,b,c), fast protons, 
deutrons and tritons (d,e,f) and fast and thermal protons (h,i,j). Each distribution was fit to three functional forms, Eq. [1] 
(solid), Eq. [T]with Co = 1 (dashed), and a second order polynomial (dotted). 



(6) are given by, 

V 

b(N gIcy ) = ^bP{N m \v)x(v,b). (4) 

v,b 

The excitation energy per residue nucleon (E*/A) was determined by the method described in [l2| by summing the 
energy measured in all charged particle fragments with a parameterized contribution from neutrons based upon . 
The thermal energy cuts used to calculate contributions to the excitation energy were the same, i.e. 30 MeV for 
protons and 49 MeV for deuterons and tritons. Therefore, E*/A is independent of the N gvcy used to determine 
centrality and for all of the systematic error evaluations except for the set of fits that specifically include the thermal 
contribution. 



IV. RESULTS 



Fig. [2] shows the excitation energy per residue nucleon as a function of the measured N gTcy (a), and the extracted 
quantities for (i/) (b) and (b) (c). The systematic error bands show the RMS variation in E*/A for the three functional 
forms, the inclusion of thermal protons and fast deuterons and tritons in the N glcy definitions, and the ±2 nib variations 
in cross-section described earlier. The systematic error bands are centered about the mean for the set of all variations, 
leading to a slight displacement from the standard analysis values that is most prominent for the pions. 
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FIG. 2: (Color online) Excitation energy E*/A vs. -/Vg rcy (a), the mean number of hadron-nucleon scatterings, {v) (b), and the 
mean impact parameter, (6) (c) plotted for three systems: p-Au (circles), 7r-Au (squares), and p-Au (triangles). Filled bands 
indicate systematic errors. 



The dependence of E*/A on N gTcy is approximately linear and nearly uniform for the three systems, with a slightly 
higher energy deposition for the antiprotons. This uniformity does not extend to the extracted mean number of 
scatterings and mean impact parameter, except in the case impact parameter dependence for the most central collisions 
((b) < 4 fm). The relation between N glcy and the calculated values of (v) and (b) for each system is a direct consequence 
of the different hadron-nucleon cross-sections. To reach the same (v), the pions must traverse a much thicker region 
of the nucleus than the protons or antiprotons. This implies a smaller impact parameter and leads to a larger number 
of fast protons emitted and therefore a greater excitation of the nucleus. This is illustrated in Fig. [3] panels (a) and 
(f) which show the fast proton multiplicity (N gToy ) vs. the extracted quantities of (v) and (6). The remaining panels 
show the dependence of other light fragments. Fig. [1] shows the dependence of the thermal particle emission on the 
extracted centrality measures. The mean fragment values in Figures [3] and |4] include efficiency corrections described 
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FIG. 3: (Color online) Mean fast particle multiplicities vs. mean number of hadron-nucleon scatterings ((^)) and mean impact 
parameter ({6}) extracted from the number of fast protons for three systems: p-Au(circle), 7r-Au(square), and p-Au(triangle) 



in 0, [13]. Note that small differences for the most central collisions are not significant relative to the systematic 
errors, which have been omitted for clarity. 



CONCLUSIONS 



We have presented the first measurements of nuclear excitation energy and light charged particle emissions as a 
function of centrality. As expected, the excitation energy increases with centrality (see Fig. 10 of [2(3]), but it is most 
directly related to the number of emitted grey particles defined as the fast protons with 30<E<350 MeV. This is 
consistent with the two stage multifragmentation process, in which a prompt intranuclear cascade depletes and excites 
the remnant nucleus which then fragments and decays. The different dependences on (v) and (b) serve to illustrate the 
importance of the collision geometry in understanding how the cascade develops from the initial primary collisions. 
This analysis will serve as an important constraint for models such as the intranuclear cascade (INC) [2l|, |22j, the 
hybrid statistical multifragmentation model (SMM) [23|, and the Boltzmann-Uehling-Uhlenbeck model (BUU) [24j . 
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FIG. 4: (Color online) Mean thermal particle multiplicities vs. mean number of hadron-nucleon scatterings ((f)) and 
mean impact parameter ((b)) extracted from the number of fast protons for three systems: p-Au(circle), 7r-Au(square), and 
p-Au(triangle) 



In particular, it will be interesting to see if the measured N grcy dependence can be reproduced, and whether the 
extracted trends in excitation energy and light fragment emission vs. (v) and (b) are borne out within the detailed 
models. 
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